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ABSTRACT 

Weak gravitational lensing by large-scale structure affects the determination 
of the cosmological deceleration parameter q Q . We find that the lensing induced 
dispersions on truly standard candles are 0.04 and 0.02 mag at redshift z — 1 
and z = 0.5, respectively, in a COBE-normalized cold dark matter universe 
with tt = 0.40, A = 0.6, H = 65km/s/Mpc and a 8 = 0.79. It is shown that 
one would observe q = —OA4to'^l and q = — 0.45lo!o3 (the errorbars are 2a 
limits) with standard candles with zero intrinsic dispersion at redshift z — 1 
and z = 0.5, respectively, compared to the truth of go — —0.40 in this case, i.e., 
a 10% error in q will be made. A standard COBE normalized fi = 1 CDM 
model would produce three times as much variance and a mixed (hot and cold) 
dark matter model would lead to an intermediate result. One unique signature 
of this dispersion effect is its non Gaussianity. Although the lensing induced 
dispersion at lower redshift is still significantly smaller than the currently best 
observed (total) dispersion of 0.12 mag in a sample of type la supernovae, 
selected with the multicolor light curve shape method, it becomes significant 
at higher redshift. We show that there is an optimal redshift, in the range 
z ~ 0.5 — 2.0 depending on the amplitude of the intrinsic dispersion of the 
standard candles, at which q Q can be most accurately determined. 

Subject headings: Cosmology: large-scale structure of Universe - cosmology: 
theory - gravitational lensing - numerical method - supernovae 
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1. Introduction 

Any interpretation of observations at cosmological scales is highly dependent on the 
density of the universe along the line of sight to the observed object. It is precisely for this 
reason that many types of observations are made, to determine the mean density along the 
line of sight. But most of the classical cosmological tests (Weinberg 1972; Peebles 1980) 
have been designed on the assumption that the global mean density p will be measured. 
At very high redshift the COBE observations tell us that the universe was very uniform. 
However, we now know that the fluctuations about this mean are large [(Sp/p) rms > 1, 
for smoothing scales less than 5/i -1 Mpc], and that even the surface mass densities along 
different lines of sight can show large fluctuations [(5E/E) rms > 1] at low redshift. Also 
gravitational instability theory tells u uss that the growth of dumpiness in the low to 
intermediate redshift range depends sensitively on some still largely uncertain cosmological 
parameters such as Qq and go- Such a non- uniform matter distribution between our local 
observing point and distant light sources will affect the properties of the sources in two 
different ways: the angular size of extended objects is changed (Gunn 1967a) and the 
apparent brightness of a source is affected (Gunn 1967b). 

The rare very large fluctuations when E « p CTit cH§ x « 10 2 E will produce gravitational 
lensing of multiple images, where the source and the observer are connected by more than 
one null geodesic, and two or more images of a background quasar (first discovered by 
Walsh, Carswell and Weyman 1979) or galaxy (for an impressive recent example see Colley, 
Turner & Tyson 1996) can be seen. This happens along a trivial fraction of lines of sight 
estimated to be a fraction ~3x 10~ 3 from observations of double quasars. Precisely the 
rarity of such events make them sensitive tests of cosmological models (|Gen et al. 1994| ; 



Kochanek 1995; Wambsganss et al. 1995| ). The much more common effect caused by 



gravitational lensing, due to much smaller fluctuations in E ("weak" gravitational lensing), 
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appears as either some shape distortion of background galaxies by large scale matter 
distribution, primarily in the outskirts of galaxy clusters (|Tyson et al. 1990| ; Miralda-Escude 



1991| ; [Kaiser fc Squires 1993| ; for a recent observational example see |Squires et al. 19961) or 



some change in the apparent brightness of background sources. In this letter we study the 
latter effect: apparent brightness changes of moderate redshift light sources caused by the 
large-scale structure, showing that weak gravitational lensing can alter the determination 
of the cosmological deceleration parameter go m significant ways. In the context of ad hoc 
assumptions for the matter distribution based plausibly on observations see e.g. Ostriker & 
Vietri (1986). Here we compute such a gravitational lensing effect in a concrete manner 
using specific models for the growth of structure. We apply it to type la supernovae as 
"standard candles" used for the determination of the cosmological deceleration parameter 
go- A complementary analytical study of this effect, based on some approximations of the 
cosmic matter distributions by simple models with adjustable parameters was recently done 
by Kantowski et al. (1995), 



2. Gravitational Lensing Magnification Distribution 

The details of our ray tracing method to study the gravitational lensing effects in both 
strong and weak regimes can be found in Wambsganss et al. (1996), and the first results 
focusing on the strong lensing effect in a standard Cold Dark Matter model have been 
published (Wambsganss et al. 1995). Here we briefly reiterate the basics. In order to study 
the lensing effects of a cosmogonic model due to the large-scale cosmic structure, large-scale 
N-body simulations are used to produce contiguous matter distributions at all redshifts. 
For the convenience of data handling (without compromising the wanted accuracy), we in 
practice approximate such a contiguous three dimensional matter distribution by a stack of 
about sixty two dimensional surface matter planes (lens planes) filling up the redshift range 
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between observer and sources at high redshift. This treatment of the matter in the universe 
in many lens planes acts effectively as three dimensional gravitational lensing. Then we 
follow a large number of light rays from an observer at z = through many lens planes 
up to a source plane at high redshift. In each lens plane the deflection of the light rays is 
determined due to all the matter in this plane. 

To quantify the effects, we use a very high resolution simulation of a cold dark matter 
model with SI = 0.40, A = 0.6, H = 65km/s/Mpc and a 8 = 0.79 (LCDM) (normalized to 
first year COBE; Smoot et al. 1992). Note that this model produces both abundances of rich 
clusters of galaxies at present epoch and large-scale structure consistent with observations 
(Bahcall & Cen 1992; Cen & Ostriker 1994; Peacock & Dodds 1994; Cen 1996; Eke, Cole & 
Frenk 1996). The simulation is run using the TPM ( |Xu 1995| ) N-body code with a box size 



of 80/i -1 Mpc. A PM grid of 512 3 points and subcell resolution extension (Tree resolution) 
of 30 are used to give a total nominal dynamic range of 15370, translating to the nominal 
resolution of 5/i _1 kpc or a true resolution of ~ 13/i~ 1 kpc. Since what we are interested 
in here is weak gravitational lensing effect, both the achieved resolution on small scales 
and the large volume of the simulation is adequate. [Note that for the present purpose 
hydrodynamic simulations (e.g. (Jen 1992 ) are not needed, for baryonic matter does not 



contribute significant amount of mass until one reaches galactic scales where cooling causes 
high condensation of the baryons.] This large-scale simulation has sufficient dynamic range 
so that we do not need a "convolution" step (Wambsganss et al. 1996). 

We fill the universe with square planes of 8/i _1 Mpcx8/i _1 Mpc comoving (with a 
thickness of comoving size 80/i -1 Mpc), produced by projecting rectangular cylinders 
of 8/i _1 Mpcx8/i _1 Mpcx80/i~ 1 Mpc, which in turn are taken from the L = 80iif _1 Mpc 
simulation boxes at different redshifts and at random angles. Each of these lens planes has 
800 x 800 pixels of size 10 x 10/i _2 kpc 2 . We follow an originally regular grid of light rays 
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through the lens planes up to various source redshifts, considering the deflection of each 
ray in each plane by all the matter of the plane. As a result we obtain the positions of 
the light rays in the source plane in a field of angular size of about (6arcmin) 2 . With the 
simple relations between these positions and those in the sky plane, we obtain properties 
like shear, effective surface mass density, and - what is important for the purposes here - 
the magnification as a function of position in the source plane (cf. [Schneider et al. 1992| ; 
Wambsganss et al. 1996| ). 

A total of 80 independent stacks of planes are sampled ("lines of sight"), each with 800 2 
positions, and we determine the magnification distribution f(fi)d/j, for sources at redshifts 
z = 0.5 and z = 1.0 (see Figure 1). The z = 1 results are shown as thick solid curves and 
the z = 0.5 results are shown as thin solid curves. The curves in the top panel describe 
the probability that a source is (de-)magnified by //, in the lower panel the cumulative 
probability is displayed: /(> //). 

Note that the median of the cumulative distribution is left of /x = 1: most sources will 
be slightly demagnified. The average magnification is of course n = 1 (i.e. the same as if all 
matter were smeared out), and a small number of sources is moderately magnified, with a 
very small high magnification tail. There is a minimum value of the magnification for each 
redshift, this corresponds basically to an "empty beam" case, in which no matter is inside 
the ray bundle. The /(> //) = (97.5%, 50%, 2.5%) points are /i = (0.951,0.983, 1.101) and 
fj, = (0.978,0.993, 1.034), respectively, at redshifts z = 1 and z = 0.5. We will use these 
magnification distributions to compute the effects on the determination of go- To set the 
context, let us now discuss the particular observations concerning the type la Supernovae. 
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3. Type la Supernovae as Standard Candles: Determination of q 

Recently, at least a dozen supernovae at redshifts of z — 0.3 or beyond have been 
detected (Perlmutter et al. 1995a,b,c; Garnavich & Challis 1996a,b), most of them of 
or consistent with type la. They can be used to determine the Hubble constant H and 
the cosmological deceleration parameter go (Perlmutter et al. 1995a, 1996; Riess, Press & 
Kirshner 1995). Currently the most distant supernova detection is at z = 0.65 (Perlmutter 
et al. 1995c). One assumes knowledge of the absolute brightness of these supernovae or 
a subset with additional characteristic features, such as "well-measured light curves" and 
"not being unusually red or spectroscopically peculiar" ([Branch fc Miller 1993j ; Vaughan 



|et al. 1995| ; Perlmutter et al. 1995a), known "minimum ejection velocities of calcium" 
(Fisher et al. 1995), or "Light Curve Shapes" (LCS; |Phillips 1993; ; |Hiess et al. 19951 ). 

For all SN la in such a subset it is assumed that they are standard candles or their 
brightness depends in a simple way on an observable parameter. The observed spread, e.g., 
of 27 "normal" SN la's as reported by Perlmutter et al. (1995a) is <jm b = 0.3 mag in B 
or om v =0.25 mag in V, with the intrinsic dispersion expected to be even smaller. Riess, 
Press, & Kirshner (1996) use a multicolor light curve shape method (MLCS) and are able 
to select a subset of type la supernovae with a dispersion of only 0.12 mag. Small intrinsic 
dispersions make SN la excellent candidates for measuring the fundamental cosmological 
parameters including Hubble constant H and the deceleration parameter q (for a recent 
review on the subject see |Branch fc Tammann 1992] ). It is plausibly assumed that there is 



no systematic trend of the intrinsic luminosity of the supernovae with increasing redshift 
(although since the mean metallicity may be expected to be redshift dependent, this 
assumption should be examined). 

Perlmutter et al. (1995a, 1996) have shown that it is possible to determine the value of 
the cosmological parameter q by detecting supernovae Type la at high redshifts (z > 0.5). 



- 8- 



We, however, will show that gravitational lensing of the standard candles by the intervening 
large-scale matter distributions causes a systematic, though not large, underestimate of the 
true go, even if the standard candles are perfect (i.e., their intrinsic luminosity dispersion is 
zero). More importantly, it places fundamental lower limits on the dispersion of the derived 

go- 

Now, assuming that type la supernovae are perfect standard candles with zero 
luminosity dispersion, we can use the magnification distribution shown in Figure 1 to 
compute the effect of gravitational lensing on making them appear to be imperfect standard 
candles. In Figure 2, the thin dotted curve indicates the relationship between redshift 
(y-axis) and distance modular, m — M (x-axis) [see equation (13.30) of Peebles 1993]. Only 
a small redshift range near z=1.0 and z=0.5, respectively, in top and bottom panels is shown 
to best read the effects. If the matter in the universe were truly smooth, perfect standard 
candles would lie on the thin dotted curve. In both panels, the three solid dots represent 
97.5%, 50%, and 2.5% of luminosity distribution (converted to distance modular) of the 
type la supernovae, when the gravitational lensing effect by large-scale matter distribution 
is included. To quantify the effect in terms of introducing uncertainties in determining g , 
we show in Figure 2 three additional curves: left thick solid curve, thin solid curve and right 
thick solid curve correspond to three different universes with q = (—0.27, —0.44, —0.49) and 
Qo — (—0.35, —0.45, —0.48), respectively, in the top and bottom panels (a flat cosmological 
model with a non-zero cosmological constant is assumed for simplicity and consistency). 
These three curves represent the best fits to the solid dots. We see that even perfect 
standard candles with zero dispersion in a 

q = -0.40 

universe will be observed/interpreted as if they reflected a 

go = -0.44tj}£ (1) 



- 9- 



universe for z s = 1.0, and a 




(2) 



universe for z s = 0.5, due to only the gravitational lensing effect by large-scale structure. 
Here the value of go is the median one, and the errorbars represent 2a limits. Note that 
the lensing induced dispersion on g is non-Gaussian, which was pointed out in an earlier 
study on the same subject by Dyer & Oattes (1988). A similar calculation by Kantowski 
et al. (1995) gave a larger, 50% effect on the median value of g , compared to a 10% effect 
found here. The difference is likely due to different models adopted here and in their study, 
in that we use an observationally favored, realistic cosmological model while their model is 
highly simplified and not observationally constrained. 



We use a spatially flat cold dark matter model with Q Q = 0.40, A = 0.6, 
H = 65km/s/Mpc and a$ — 0.79, to compute the effect of gravitational lensing by 
large-scale matter distribution on the determination of go- This model agrees with 
observations on the abundances of rich clusters of galaxies and large-scale structure as 
well as is normalized to COBE. It is shown that, even with perfect standard candles with 
zero dispersion, the observed dispersion, after being convolved with the lensing effect, is 
sizable with a dispersion of 0.04 and 0.02 mag at redshifts 1 and 0.5, respectively, in the 
LCDM universe. Moreover, since the distribution is asymmetric, the median observed 
value of g is lowered with resulting values of g = — 0. 44lo!o5 an d Qo = — 0. 451q!o3 f° r 
z s = 1.0 and z s = 0.5, respectively (the errorbars are 2a limits), compared to the truth of 
g = —0.40 in this case, i.e., a 10% error in g will be made. The non-Gaussian signature 
clearly distinguishes this lensing dispersion from other possible broadening of the supernova 
la luminosity distribution. It is worth noting that, although these dispersions are still 



4. 



Discussion 
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smaller than the current observed total dispersions of some carefully selected subsamples 
of type la supernovae at lower redshift [for example, 0.12 mag from Riess et al. (1996) for 
SN la at redshift between 0.007 < z < 0.1, 0.21 mag from Perlmutter 1996| for SN la at 



redshift between 0.35 < z < 0.5], it will become increasingly important for supernovae (or 
any other standard candles with comparable or small intrinsic luminosity dispersions) at 
higher redshift. Similarly, the lensing induced dispersion will become more important if the 
observational errors can be further reduced. 

The size of this dispersion, especially the upward dispersion of go, will depend on the 
degree of dumpiness of the large-scale matter distribution as well as on the values of Qq 
and go- We also examined the standard Qq — 1 cold dark matter (SCDM) normalized to 
the first year COBE (cr 8 = 1.05), and found that the lensing induced dispersions are three 
times as large as those in the LCDM model presented here. This is not surprising since it is 
known that the SCDM model overproduces the present day rich cluster numbers by a large 
factor (|Bahcall fc Cen 1992| ; |White, Efstathiou, fc Frenk 1993|) , while the adopted LCDM 



model appears to agree with most observations QBahcall fc Cen 1992| ; Ostriker & Steinhardt 



1995; |Cen 199q ; Eke et al. 1996). We estimate from other calculations that a hot and cold 
dark matter model (HCDM) with Qq = 1 would produce an intermediate result, perhaps 
twice as large as the a = 0.04 quoted here for the LCDM model. 

Even though the differences in the geometry of the universe for cosmological models 
with different go are increasingly larger at higher redshifts (precisely for this reason it 
is generally thought that higher redshift sources are better suited to determine g ), the 
distortions/dispersions caused by gravitational lensing on go also become larger at higher 
redshifts. Furthermore, gravitational lensing due to matter fluctuations on scales smaller 
than what is considered here, such as lensing by individual galaxies and microlensing by 
stars or brown dwarfs (Linder, Schneider, & Wagoner 1988; Rauch 1991), will further 
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increase the amplitudes of the dispersion. Therefore, it is not clear whether one gains or 
loses by resorting to very high redshift sources, or what the "optimal" redshift range for 
such studies is. 

We show, in Figure 3, rj(z) = d{ *° M) ^ AMf nstrinsic + Am^(z) as a function of 
redshift z, where AM int rinsic denotes the intrinsic magnitude dispersion of the standard 
candles, and Ami ensing (z) represents the lensing induced dispersion in the apparent 
magnitude (we have assumed that lensing induced magnitude dispersions due to other, 
smaller scale structures will increase those due to large-scale structure alone by a factor 
of two). T)(z) quantifies the dispersion in the determined q given the observed dispersion 
in the apparent magnitude. Four cases for AM intrinsic = (0.20,0.12,0.06,0.03) mag are 
plotted (again, a flat cosmological model with a non-zero cosmological constant is assumed 
for simplicity and consistency). We see that one benefits greatly by going from redshift 
zero to about redshift 0.5; going further in redshift does not grant a significant gain. In 
fact, we find that the optimal redshift of SNe la for determining q is (2.0, 1.4,0.70,0.45), 
respectively, for the four assumed AM intrinsic values. Note that this effects applies to any 
set of standard candles; it is not limited to Type la supernovae. 

Finally, we point out that, if there exists a set of truly standard candles with a 
negligibly small intrinsic luminosity dispersion, it should be possible to use them to directly 
probe the gravitational lensing effect discussed here. But, on the other hand, we warn that 
even quite tiny systematic changes in either the intrinsic luminosity of the sources with 
increasing look-back time or a redshift dependent variance could easily give quite incorrect 
results for both q and the detected lensing. 

The work is supported in part by grants NAG5-2759, AST91-08103 and ASC93-18185. 
Discussions with P. Hoflich, C. Hogan, B. Kirshner, G. Lake, S. Perlmutter and C. Stubbs 
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Figure Captions 

Fig. 1.— Magnification distribution due to gravitational lensing by large-scale matter 
distributions in a A-dominated, COBE (first-year) normalized flat cold dark matter model 
for sources at redshifts z = 0.5 (thick lines) and z = 1.0 (thin lines). The curves in the top 
panel describe the probability /(//)that a source is (de-) magnified by /i, in the lower panel 
the cumulative probability is displayed: /(> ji)dji. 

Fig. 2.— The thin dotted curve shows the location of perfect standard candles in the 
distance modulus (m — M) - redshift plane for a go — —0.40 universe. The three solid dots 
represent 2.5%, 50%, and 97.5% of luminosity distribution (converted to distance modular) 
of the type la supernovae at z — 1 (top panel) and z = 0.5 (bottom panel), respectively, 
after the large-scale structure lensing effect is considered. The left thick solid curve, thin 
dotted curve and right thick solid curve correspond to three different universes with q' s as 
indicated, (a flat cosmological model with a non-zero cosmological constant is assumed for 
simplicity and consistency). These three curves represent the best fits to the solid dots. 

Fig. 3. — We show rj(z) (see text), which is proportional to the dispersion in the determined 
g , as a function of redshift z. Four cases for the intrinsic magnitude dispersion of the 
standard candles AM int rinsic = (0.20,0.12,0.06,0.03) mag are plotted (a flat cosmological 
model with a non-zero cosmological constant is assumed for simplicity and consistency). 
Note that, for a given intrinsic magnitude dispersion of a sample of standard candles, the 
minimum of r)(z) denotes the optimal redshift at which go can be most accurately determined, 
using any type of standard candles. 
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